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SUM4AKI 


* 


An  investigation  is  undertaken  to  (l)  reevaluate  the  elastic  theory 
results  obtained  by  previous  investigators  for  both  perfect  and  initially 
imperfect,  axially  compressed  circular  cylindrical  shells,  (2)  develop 
new  elastic  stability  results  describing  the  postbuckling  behavior  of 
curved  plates,  and  (3)  present  criteria  based,  in  part,  upon  the  curved 
plate  results  which  can  be  used  to  size  and  locate  stringers  and  rings 
such  that  the  initial  buckling  load  will  occur  in  the  neighborhood  of 
the  classical  buckling  load  for  a  cylinder  representing  a  minimlzed- 
weight  design. 
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FOREWORD 


The  vork  reported  herein  constitutes  a  portion  of  a  continuing  effort 
being  undertaken  at  Stanford  University  for  the  U.  S.  Army  Aviation 
Materiel  Laboratories  under  Contract  DAAJ02-68-C-0035  (Task 
1F162204A17002 )  to  establish  accurate  theoretical  prediction  capability 
for  the  static  and  dynamic  behavior  of  aircraft  structural  components 
utilizing  both  conventional  and  unconventional  materials.  Predecessor 
contracts  supported  investigations  vhich  led,  in  part,  to  the  results 
presented  in  references  1 6,  19,  24,  and  25* 
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approaches  zero.  Meanwhile,  by  utilizing  a  more  exact  set  erf  strain- 
displacement  relations  than  von  Kanaan-Donnell,  Mayers  and  Hehfield,  in 
1965/  bad  obtained  radius-to-thickness -dependent,  load-end  shortening 
curves  for  the  Kempner  deflected  shape  and  showed  that  the  von  Karman- 
Donnell -theory  solution  obtained  by  Kempner  is  reached  only  as  the 
thickness-to-radius  ratio  of  the  cylinder  approaches  zero.  As  is  demon¬ 
strated  in  reference  15,  the  deflected  shape  of  the  converged  potential 

energy  solution  of  the  von  Karman-Donnell  equations  approaches  the 

17 

developable-surface  shape  derived  by  Yoshiiura.  Thus,  in  view  of  the 
findings  of  Hoff  et  al.1^  and  Mayers  and  Rehf ield,1^  it  became  evident 
that  elastic  analysis  postbuckling  curves  for  perfect  shells  were  of  no 
practical  significance.  Attention  was  then  turned  to  the  elastic 
behavior  of  imperfect,  axially  compressed,  circular  cylindrical  shells  in 
attempts  to  more  accurately  predict  the  initial  buckling  or  maximum  load. 
The  major  effort  in  this  area  is  included,  also,  in  the  work  of  Madsen 
and  Hoff  reported  in  reference  12. 

However,  deapite  the  extensive  and  important  information  concerning  the 
elastic  stability  and  postbuckling  behavior  of  axially  compressed  cyl¬ 
inders  included  in  references  1-17,  there  is  still  no  theoretical 
analysis  which  will  predict  the  initial  elastic  buckling  load  of  a 
practically  fabricated,  thin,  axially  c expressed  cylinder.  More  impor¬ 
tant,  since  pure  monocoque  shells  are  the  exception  rather  than  the  rule 
in  shells  designed  for  compression  loading,  there  is  no  theoretical 
analysis  which  can  (1)  predict  the  initial  buckling  load  of  a  thin 
stiffened  cylinder  even  if  the  behavior  is  elastic  and  (2)  provide  design 
criteria  whereby  the  stiffened  cylinder  will  reflect  an  optimized  design. 

It  is  the  purpose  of  the  present  investigation  to  (l)  reevaluate  the 
elastic  theory  results  obtained  by  previous  investigators  for  both 
perfect  and  initially  imperfect,  axially  compressed  circular  cylindrical 
shells,  (2)  develop  new  elastic  stability  results  describing  the  post¬ 
buckling  behavior  of  curved  plates, and  (5)  present  criteria  based  in 
part  upon  the  curved  plate  results  which  can  size  and  locate  stringers 
and  rings  such  that  the  initial  buckling  load  will  occur  in  the 


2 


J 


4 

/ 


neighborhood  of  the  classical  buckling  load  for  a  cylinder  representing 
a  minimized -weight  design.  All  three  objectives  have  been  met;  the 
findings  are  presented  in  the  form  of  equations  and  graphs.  The  results 
of  the  analyses  and  the  computations  are  presented  in  the  body  of  the 
report,  and  the  details  of  the  analyses  and  computational  procedures  are 
given  in  a  series  of  appendixes. 
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GENERAL  TBEdDf 


STRAIS-DISHiACaCBIT  RELATICRS 

The  strain-displacement  relations  used  herein  are  the  small  strain, 
large  displacements -moderate  rotations  set  derived  by  Doncell  vith 
modifications  to  account  for  initial  deviations  from  uniform  circularity 
along  the  shell  length,  vith  denoting  the  initial  imperfection 

shape,  v  denoting  the  additional  deflection  due  to  bending  of  the  shell 


during  lead  application,  and  v  denoting  the  sum  of  v  and  v; 

X  v 

is. 


that 

(1) 


The  strain-displacement  equations  for  the  middle  surface  of  the  shell 
(see  Figure  l)  are 

0  A  1  ,  2  2  v  \ 

€x  =  u'x  +  2  (wT.  ~V0,  >  I 

*x  *x  I 

€°  -  V  +  i  (V2  -,2  )  -  (VV0>  1(2 

y  V,y  2  'WT,  V0,  '  R  / ' 

y  y  I 

Vs  u’y  +  V'x  +  vT,xwT,y  -  vo,xvo,y  ) 

In  view  of  equation  (1),  equations  (2)  are  rewritten  as 


0  12 
ex  =  u'x  +  2  w'x  +  w,xv0. 


0  12  w 

y  ’y  2  *y  ’y  0,  R 

0 

7  =  u,  +  v,  +  w,  w,  +  vn  w,  +  w,  v_ 

'xy  ’y  ’x  ’x  ’y  0,x  fy  ’x  0,y 


The  bending  strains  are 


% 


POTENTIAL  ENERGY 


The  total  potential  energy  stored  in  the  elastic  shell  under  prescribed 
unit  end  shortening  e  (see  Figure  1)  is  composed  of  the  extensional  and 
bending  strain  energies;  it  is  given  by 


U  = 


E 

2(1- v2) 


j.L/2  2rtRr  t/2 
J-L/2^0  "-t/2 


dxdydz 


(8) 


Substitution  of  equations  (5)  in  equation  (8)  and  subsequent  integration 
over  the  constant  thickness  t  of  the  shell  vail  yields 


o  =  f  L/2f  4»? «„  +v,  4.? v.  - 
2d-2)  1  x  2  x  ’x  °'x  y 2  'y  y  \  R 

-  Sd-vJKu,/!*?^,^  >(v,y4*?yw,y«0,  -  |)  - 


+  '-/'A  v’y+u’xvO,  > 
x  y 

r\  i»  h/2  2ltR  p  p 

+  2  J-L/s-o  -2(1-v)(“’xxv'yy-u'xy))dxdy 

vhere  D  =  EIt^/l2(l-v2) . 


(9) 


In  accordance  vith  the  potential  energy  principle,  U  must  have  a 
stationary  value  vith  respect  to  admissible  variation  in  the  degrees  of 
freedom  u,  v,  and  v  in  order  for  equilibrium  to  be  achieved.  As 
derived  in  Appendix  I,  the  variation  leads  to  the  equilibrium  equations 


N 


+  N 


=  0 


N 


+  N.. 


=  0 


N 

W  v  -  X  -  ,J  -  fy’yy-'o,^ 

-  '  0 


(10) 
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and  the  accompanying  boundary  conditions  at 


=0  or  5v  =  0  (13) 

ALTERNATE  FORM  OF  TOTAL  POTENTIAL  ENERGY  FUNCTIONAL 

With  the  use  of  the  middle  surface  strain-displacement  equations  (3),  the 
total  potential  energy,  equation  (9),  can  be  expressed  in  the  form 
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0  =  -JV  r  1/8  f  8ltR[(c0+€0)^-2(l-v)(e0€°-i  7°2)]dxdy 
2(l-v")  -L/2"0  x  y 

D  f  ^  r  2  2 

+  1  -Ll/2-0  -2(1-v,("'«v'yy-v'xy))'ix^ 


Next,  vith  the  constitutive  relations  (7)*  written  for  the  middle  sur¬ 
face,  the  strains  e°,  e°,  and  7°  can  be  replaced  to  give  the  total 

x  y  xy 

potential  energy  as 

u  -  k  J  Jj0  <<"K,2'2(1+v,(^-4))toiy 

D  P  p  2tcR  _  _ 

+  5j.L/J0  -2(1-v)<w'«','yy-v'xy)]ax^  <15> 

The  first  two  equilibrium  equations  (10)  car  be  satisfied  identically  by 
the  Airy  stress  function  F,  defined  as 


N°/t  =  a°  =  F. 


N°/t  =  <j°  =  F, 

Y  y  XX 

N°  /t=  T°  =  -F, 
xy'  xy  'xy 


finally,  the  introduction  of  the  relationships  between  the  middle -surface 

stresses  cr°,  ct°,  and  t°  and  the  Airy  stress  function  F  defined  by 
x  y  xy 

equations  (16)  permits  the  elimination  of  the  stress  components  from 
equation  (15);  thus  the  desired  alternate  form  of  the  total  potential 


energy  is 


.  a  L/2  ^2*®  p  p 

u  =  ~  [(F,  +F,  )  -2(l+v)(F,  F,  -F,  )]dxdy 

2Ej_Ly2JQ  LV  'yy  'xx'  v  /v  'yy  'xx  'xy7 

n  (*  ^/2  f*  ?  2 

+  2  J_l/2*'o  -2(1-'')('''xx''’yy-’,'xy)]dx,iy  (17) 


The  total  potential  energy-  consists  now  of  only  two  unknown  quantities: 
the  stress  function  F  and  the  additional  radial  displacement  v  of 
the  shell.  Vith  middle-surface  equilibrium  automatically  satisfied 
through  introduction  of  the  stress  function,  the  variational  principle 
need  be  applied  only  vith  respect  to  v  to  achieve  equilibrium  out  of 
the  middle  surface.  However,  it  is  necessary  first  to  establish  a 
relationship  between  F  and  v  in  order  to  present  the  total  potential 
energy  in  terms  of  v  alone.  From  displacement  compatibility  in  the 
classical  theory  of  elasticity,  it  is  known  that  the  middle  surface  dis¬ 
placements  u  and  v  must  satisfy  the  condition 

(u’x>’yy  +  '  (u'y+T’x)’xy  =  0  (lB) 

With  the  strain-displacement  equations  (3)  and  the  constitutive  relations 
(7),  written  for  the  middle  surface,  the  stress-displacement  compatibil¬ 
ity  equation  becomes 

0  0  0 
cr  or  T 

V  E;'yy  v  E ’’xx.  K  E  "xy 


+  2wrt  w,  -v_  w,  -v„  w, 

0,  'xy  0,  'yy  0,  'xx 
’xy  ^  ’xx  ’yy 


Then,  by  introducing  the  relations  between  the  middle-surface  stress 
components  and  the  Airy  stress  function  F  of  equations  (16),  the  com¬ 
patibility  equation  of  von  Khrman -Donnell  shell  theory  is  obtained  for 
linear  elastic  behavior  as 


k  i  2  l 

V  F/E  =  v-  -  w,  w,  -  —  w,  +  2vn  w, 

'  xy  'xx  'yy  R  'xx  0,  ': 


w,  -  v~  w, 
°'xx  °’yy  '« 


For  any  specific  case,  in  order  to  eliminate  the  stress  function  F, 
the  procedure  is  to  solve  equation  (20)  for  F  in  terms  of  w.  Then, 
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vith  F  known,  the  total  potential  energy  becomes  a  function  of  v 
alone.  This  procedure  is  not  a  simple  one,  since  equation  (20)  lends 
itself  to  ready  solution  only  when,  for  example,  trigonometric  series 
are  used  to  express  v.  Further,  it  should  be  pointed  out  that  the 
procedure  is  highly  restricted  1  general  to  materials  vith  constant 
elastic  moduli  (independent  of  the  coordinates  of  the  body)  in  view  of 
the  operations  implicit  in  deriving  equation  (20)  from  equation  (19)* 


METHOD  OF  SQLUTIC1I 


With  the  total  potential  energy  applicable  to  both  imperfect  a*t«ny 
compressed  curved  plates  and  shells  given  by  equation  (17)  in  terns  of 
the  deflection  w  and  the  stress  function  F,  the  direct  method  of  the 
calculus  of  vs^iations  ( ftayleigh-Ritz  method)  is  utilized  to  establish 
equilibrium- state,  average  stress-end  shortening  curves  for  various 
imperfection  patterns. 

CIRCULAR  CYLINDRICAL  SHELL 

To  study  the  various  facets  of  the  compressed  circular  cylindrical- shell 
problem  and  to  establish  a  better  understanding  of  solutions  obtained 
through  von  Karman-Donnell  theory,  in  conjunction  with  the  potential 
energy  method  and  trigonometric  series  to  represent  the  radial  displace¬ 
ment  patterns,  both  initial  and  additional,  several  cases  are  investi¬ 
gated. 

General.  Case  of  Present  Analysis 

A  reasonable  imperfection  pattern  for  a  long  circular  cylindrical  shell 
is 


12.  A° 
t  11 


cos 


ra  ity 
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+  A, 
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cos 


2jix 


(21) 


This  expression  represents  a  modification  of  the  chessboard  pattern  into 
which  s  perfect  cylinder  buckles  with  infinitesimal  anplitudes  according 
to  the  classical  theory  (see  references  1,  2,  and  3).  It  is  also  both 
the  experimentally  observed  ouckle  pattern  of  a  compressed  shell  (see 
reference  18)  and  the  modai  pattern  of  a  long  thin  cylinder  undergoing 
free  vibration  with  finite  amplitudes  (see  reference  19).  Thus,  the 
initial  pattern  selected  fcr  w^  is  obviously  one  to  which  a  compressed 
shell  is  particularly  sensitive. 

To  provide  for  the  growth  of  the  initial  pattern  wQ  upon  application 
jf  end  shortening,  the  expression 
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(22) 


cos 


cos 


2*x 


Is  assumed.  Equation  (22)  ensures  that  the  imperfection  pattern  w^ 
retains  its  original  wavelengths  hut  is  free  to  grow  in  amplitude 
through  A*^  and  A|q. 

The  huckle  mode  for  the  shell  is  taken  to  he 


w  * 
t  "  *11 


cos 


itx  ny  2  4  2sx  ,  . 

r  COS  X  +  *20  C0S  X”  +  * 

x  y  x 


02  cos  T1  +  *oo  +  T 
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(23) 


It  is  to  he  noted  that  except  for  A__,  \  ,  and  \  ,  all  of  the 

Xq  y0 

parameters  in  equation  (23)  ere  arbitrary;  they  are  to  be  determined  by 
application  of  the  variational  principle.  That  is,  a  stationary  value 
of  the  total  potential  energy,  equation  (17),  is  sought  with  respect  to 
AJj,  A5q,  Ajj,  A 2qj  A Q2,  Xx,  and  Xy.  Since  Aqq  actually  vanishes 
before  variation  of  the  total  potential  energy,  it  is  determined 
Independently  from  enforcement  of  the  periodicity  condition  governing 
the  circumferential  displacement  v. 

Of  course,  before  the  variation  is  carried  out,  equation  (21)  and 
equation  (23)  are  substituted  into  equation  (20)  to  determine  F.  Such 
a  procedure  leads  to 


F 
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-  §r +  -s— +  v 


'XX 


xy  'xx  fyy  R 


0,  ’xy 
’ xy 


-  w„  w,  -w_  w,  ) 
0,  *yy  0,  ’xx’ 
’xx  ’yy 


(2k) 


where  cr  has  been  introduced  to  provide  for  an  average  compressive 
stress  acting  on  each  shell  cross  section.  The  actual  expression  for  F 
obtained  by  the  operations  indicated  in  equation  (24)  is  illustrated  in 
Appendix  II.  The  average  compressive  stress  is  developed  also  in 
Appendix  II  for  the  same  case  and  is  shown  to  be 
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cr  t,  2_,^11  .  .2  x  .  2_  A  .  .,2  .  *11*11  .  ..  x, 

E  =  6  “  rt-BT  +  *20}  +  Wr  +  *&  +  ~ir”  +  2*20*£o);l 


(25) 


where  e  is  the  prescribed  unit  end  shortening  and  thus  Is  the  forcing 
function  for  the  problem. 


A  typical  average  stress-unit  end  shortening  curve  for  an  assumed 


inperfection  pattern  wherein  A^  = 
and  Tj  -  s^Rt/X2  =  0.25  is  given 

0  Yq 

for  pn  and  rjn  correspond  to  the 
u  ^  13 

portion  of  the  Kempner  curve  when 

„  .  8  , 

von  Xarman-Tsien  Solution 


1,  4)  a  °-25,  Ko  -  *y  Ax  =  0.52, 
in  Figure  2.  The  values  selected 
p  and  tj  values  of  the  stable 
eR/t  =  0.70. 


The  first  solution  of  the  von  KArman -Donnell  theory  for  perfect  shells 
appears  in  reference  8.  The  main  contribution  of  this  work  lies  in  the 
qualitative  rather  than  quantitative  results,  since  the  authors  show 
the  existence  of  equilibrium  states  with  large  as  well  as 
(infinitesimal)  displacements.  The  former  imply  that  equilibrium  stress 
states  exist  far  below  the  cne  corresponding  to  classical  buckling.  This 
important  finding  became  the  plausible  explanation  for  axially  compressed 
thin  cylinders  in  actual  test  to  achieve  only  a  small  fraction  of  the 
classical  buckling  load. 


The  constrained  deflection  utilized  by  von  Kannan  and  Tsien  is  equivalent 
to  the  deflection  function  of  equation  (25)  with  w^  =  0  and  A^  =  A^q- 
Typical  von  Ka’rman-Tsien  average  stress-end  shortening  curves  are  plotted 
in  Figure  5  for  p  =  X ^/Xy  =  ^  ^  several  values  of  rj  =  j^Rt/X2 
(0.676,  0.400,  0.255,  and  0.169).  A  curve  obtained  from  the  present 
analysis  is  plotted  also  in  Figure  5  for  p  =  X ^/X^  =  1,  A^  =  A^0, 
and  tj  absolutely  free. 

Kempner^  Solution 

The  first  unconstrained  solution  of  the  von  Kannan-  Donnell  theory  for 
perfect  axially  compressed  thin  shells  is  given  in  reference  15.  This 
solution,  for  the  deflection  function  assumed  (equation  (25)  with 
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v*  s  o  and  Ajj,  Ag^,  A^,  n,  and  ij  completely  arbitrary)  represents 

a  true  equilibrium  solution,  since  all  radial  deflection  parameters  are 

peraitved  to  vary  in  the  application  of  the  variational  principle.  The 

calculated  points  obtained  by  Keapner^  to  construct  his  average  stress- 

,  * 

end  shortening  curve  are  shown  in  Figure  4. 

Madsen-  Hoff12  Solution 

In  references  4  and  12,  Hoff  et  al.  point  out  the  basic  error  of  Donnell 
and  Van  appearing  in  reference  U,  that  of  actually  permitting  the 
imperfection  pattern  amplitude  and  wave  shape  to  vary.  In  reference  12, 
the  authors  avoid  the  error  by  assuming  an  imperfection  shape  and  per¬ 
mitting  only  the  growth  of  the  imperfection  amplitude  to  vary.  Their 
assumed  initial  imperfection  shape  and  arbitrary  buckle  pattern  are 
equivalent  to  equation  (21)  and  equation  (23)  of  the  pres  nt  work  with 
A^  «  0,  Ag£  **  0,  and  A^g  «  0.  The  average  stress-end  shortening 
curves  obtained  on  this  basis  are  shown  in  Figure  3  for  A^,  =  0,  0.1, 
and  0.17.  Also  shown  in  Figure  5  are  the  results  of  the  present  study 
for  Aqj,  f  0,  ly  =  Xy  ,  A^  =  0.1,  and  0.17. 

CTROILAHMr  CURVED  PLATE 

Most  of  the  research  that  has  been  conducted  in  connection  with  under¬ 
standing  the  so-called  "perplexing  behavior"  (see  reference  5)  of  thin 
shells  in  conpreseioa  has  had,  in  reality,  the  goal  of  providing  criteria 
for  conventional  shell  structures  used  in  aerospace  design  to  withstand 
direct  compression  and  bending  loads.  Such  structures  appear,  in 
general,  in  stiffened  rather  than  unstiffened  form;  hence,  theoretical 
studies  on  pure  monocoque  shells  under  direct  compression  and  bending 
loads  have  only  limited  value  from  the  practical  viewpoint. 

A  stiffened  shell  of  the  noncomposlte  type  may  be  viewed  as  being  made 
up  of  a  series  of  curved  plates  separated  by  stringers  and  rings.  Thus, 
it  is  of  some  importance  to  study  the  behavior  of  curved  plates  as 


Professor  Joseph  Kempner  kindly  provided  the  calculated  points  from 
which  the  continuous  curve  in  reference  13  was  plotted. 
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discrete  structural  elements  and  also  to  establish  a  knowledge  of  the 
transition  from  the  well-understood  behavior  of  flat  plates  to  the  not- 
so-vell-understood  behavior  of  complete  cylinders  as  the  curvature  is 
gradually  increased. 

The  shell  theory  presented  in  GENERAL  THEORY  and  Appendix  I  is  equally 
applicable  to  curved  and  flat  (R  -*  •)  plates  singly  by  changing  the 
limits  of  integration  in  the  circumferential  direction  in  the  potential 
energy  expression  (equation  (17))  and  the  boundary  conditions  (equations 
(ll)-(l3))  from  y  =  0,  2itR  to  y  =  -  77.  naturally,  the  same 

changes  would  occur  in  the  developments  of  Appendix  I.  A  long  curved 
plate  element  is  shown  in  Figure  1,  along  with  the  sign  convention  used 
in  the  analysis. 

For  a  first  approximation  to  the  behavior  of  a  circularly  curved  plate 
under  prescribed  end  shortening,  the  radial,  deflection  function  is  taken 
as 


^  =  A^  cos  ^  cos  ^  +  A„  cos  ^  ~  =  1,3,5...  (26) 

x  y  y 

The  first  term  is  sufficient  to  establish  quite  accurately  the  average 
stress-end  shortening  curve  for  a  long  rectangular  flat  plate  (R  -*  •) 
with  pinned  edges  until  the  end  shortening  exceeds  the  critical  end 
shortening  by  an  order  of  magnitude  (see  reference  20).  For  shallow- 
curved  plates  (Z  =  b'/Rt  <  10),  the  second  term  in  equation  (26)  is 
sufficient  to  provide  reasonably  accurate  average  stress-end  shortening 
curves  dependent  upon  Z.  Through  the  use  of  equation  (26)  and  the 
total  potential  energy  given  by  equation  (17)  (with  suitable  integration 
limits  for  the  circumferential  direction),  application  of  the  variational 
principle  leads  to  the  average  stress -end  shortening  curves  plotted  in 
Figure  6,  when  p  *  \  f\  =1,  N  =  b/l  =  1,  and  Z  =  10  and  20, 
respectively.  Also,  in  Figure  6,  are  shown  corresponding  results  for 
an  improved  approximation  to  the  radial  deflection  function;  that  is. 
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I =  C06  r cos  ? +  Aoi cos  ? +  *21 cos  x2  C06  ?  (27) 

x  y  X  y 

The  details  of  the  curved  plate  analysis  are  given  In  Appendix  m. 

It  is  to  be  noted  that  the  behavior  of  a  shallow  curved  plate  (Z  <  10) 
is  slailar  to  that  of  a  flat  plate;  thus,  the  curved  plate  is  relatively 
insensitive  to  Initial  deviations  (see  reference  21). 


RESULTS  AKD  DISCU3SI0B 

CIRCULAR  CnJKDRICAL  SHELL 

An  investigation  of  the  bending,  buckling,  and  postbuckllng  of  initially 
inperfect  circular  cylindrical  shells  under  axial  compression  has  been 
undertaken  to  establish  a  clearer  understanding  of  the  many  solutions 
obtained  through  the  use  of  von  Karman-Donnell  shell  theory,  the 
potential  energy  method,  and  trigonometric  series  to  represent  the 
radial  displacement  patterns. 

General  Case  of  Present  Analysis 

The  most  general  solution  of  the  present  study  is  based  upon  the  assumed 
initial  and  additional  radial  displacement  functions  given  by  equations 
(21)  and  (23).  The  assumed  waveform  differs  from  that  of  Kempner  only 
in  the  stipulation  of  an  imperfection  pattern  and  the  allowance  for  its 
growth  in  amplitude.  With  wQ  and  w*  set  equal  to  zero,  the  stable 
branch  of  the  average  stress-end  shortening  curve  first  obtained  by 
Kempner  with  5  free  parameters  is  duplicated  in  Figure  2  by  the  lower 
curve.  The  upper  curve  represents  a  7 -free -parameter  solution  obtained 

from  the  present  analysis  with  A?,  =  1,  A?n  =  0.25,  un  =  X  /X  =  0.52, 

2/2  _u'  u  ^0  XQ 

and  t{q  =  it  Rt/X^  =  0.25.  The  values  of  and  tj0  correspond  to 

those  of  Kempner  when  eR/t  =  0.70.  The  curve  is  of  interest  because  it 

shows  that  unlike  average  stress-end  shortening  curves  for  postbuckled 

n  22 

shells  as  presented  by,  for  example,  Donnell  and  Wan,  Loo,  and 
12 

Madsen  and  Hoff,  any  perfect  body  solution  is  a  lower  rather  than 
upper  bound  to  the  inperfect  body  behavior  when  the  von  Kaiman-Donnell 
theory  is  utilized  in  combination  with  trigonometric  series  and  the 
Rayleigh-Ritz  procedure.  A  more  meaningful  interpretation  of  the  two 
curves  in  Figure  2  is  presented  later  in  connection  with  the  detailed 
discussion  of  the  Kempner1^  solution  as  reevaluated  in  the  present 
analysis. 

von  Karman-Tslen  Solution 

Although  the  von  Karman-Tsien  solution  was  undertaken  in  a  qualitative 
fashion  only,  mainly  to  offer  a  plausible  explanation  for  the  failure  of 
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compressed  shells  In  experiments  to  reach  more  than  a  sue  11  fraction  of 
the  classical  buckling  load,  many  subsequent  investigators  have  attrib¬ 
uted  quantitative  results  to  von  Karaan-Tsien  which  the  authors  them¬ 
selves  nsver  claimed. 

Thus,  as  a  first  step  in  evaluating  the  existing  work  on  the  post- 
buckling  behavior  of  imperfect  shells,  the  von  Karaan-Tsien  solution  was 
reviewed  in  the  light  of  the  present  analysis  approach.  On  the  basis  of 
equation  (23),  with  w*  =  0  and  =  k^Q>  -be  curves  plotted  in 
Figure  3  were  obtained.  The  solid  curves  Id  the  postbuckling  region 
represent  identically  the  results  obtained  by  von  Karmen  and  Tsien  for 
(i  -  1  and  the  designated  values  of  tj.  The  dashed  curve  is  the  result 
of  the  present  analysis  when  p  =  1  and  tj  is  a  free  variable.  Of 
importance  is  the  fact  that  for  p  =  1,  the  lowest  postbuckling  stress 
*  a  wrlnltmim  potential  energy  basis  is  given  by  oB/Et  =  0.243.  The 
.lowest  postbuckling  stress  obtained  by  von  Karmen  end  Tsien  appears  on 
the  r)  =  0.255  curve  and  is  given  by  crK/Efc  *  0.194.  Leggett  and  Jones'^ 
noted,  with  respect  to  the  results  of  von  Karmen  and  Tsien,  that  the 
lowest  postbuckling  stress  obtained  for  a  fixed  p  and  discrete  values 
of  rj  is  not  necessarily  that  corresponding  to  a  minimum  energy  at  the 
fixed  value  of  p.  This  can  be  illustrated  in  Figure  3  by  noting  that 
only  certain  points  on  each  tj -curve  correspond  to  points  on  the  true 
minimum  energy  curve  for  the  deflection  assumed  and  p  =  1.  The 
important  fact  is  that  the  minimum  stress  obtained  from  the  envelope  of 
the  discrete  rj-value  curves  is  lower  than  the  postbuckling  stress 
corresponding  to  the  minimum-energy  solution  for  the  fixed  value  of  p. 
In  other  words,  in  systems  which  possess  unstable  branches  of  the  post- 
buckling  curves,  it  is  not  permissible  to  obtain  the  lowest  stress  by 
minimizing  the  load  with  respect  to  the  buckling  parameter?.  It  is 
interesting,  nevertheless,  that  with  the  same  deflection  function 
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and  A20  =  A^,  the  von  Karman-Tsien  minimum  postbuckling  stress  of 
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aR/Et  =  0.194  (obtained  for  n  =  1  and  tj  =  0.255)  compares  almost 

Identically  vith  the  minima  postbuck  ling  stress,  crp/Et  =  0.195, 

9  10 

obtained  by  Leggett  and  Jones  and  Michielsen  on  the  basis  of  free 
variation  vith  respect  to  both  p  and  tj. 

Keener  Solution 

Although  a  portion  of  the  stable  branch  average  stress-end  shortening 
curve  of  the  Kempner  solution  has  been  shown  in  Figure  2,  the  solid- 
curve  result  presented  by  Kempner  in  reference  13  is  shown  in  Figure  4 
in  terms  of  the  actual  calculated  points  (see  footnote  on  page  14).  The 
solid  curve  of  Figure  4  is  the  result  obtained  in  the  present  analysis 
with  the  aid  of  a  high-speed  computer.  The  desk-calculator  results 
obtained  by  Kempner  in  a  remarkable  effort  are  quite  accurate  along  the 
stable  branch  of  the  average  stress-end  shortening  curve.  However,  even 
though  the  remaining  points  appear  to  indicate  that  the  postbuckling 
curve  closes  to  the  bifurcation  point,  the  very  accurate  numerical 
results  of  the  present  analysis  show  that  (l)  the  left -hand  portion  of 
the  lower  branch  of  the  solution  cannot  be  extended  economically  for 
afi/Et  >  0-334,  a  Burroughs  B-5500  computer  notwithstanding,  and  (2) 
another  solution  of  completely  unstable  character  exists  for 

ffR/Et  >  0.387. 

The  loop  appearing  below  the  bifurcation  point  is  not  a  new  result. 
Thielemann,  in  Figure  10  of  reference  23,  notes  that  Donnell's  solution 
yields  the  same  phenomenon.  However,  in  the  Donnell  and  War^’  study 
referred  to  by  Thielemann,  the  authors  mention  only  that  success  was 
achieved  in  finding  a  solution  emanating  from  the  bifurcation  point  down¬ 
ward  to  ofi/Et  =  0.45.  The  unstable  solution  moving  upward  from  the 
bifurcation  point  is  considered  a  new  result  in  view  of  its  absence  from 
the  literature. 

Another  new  result  of  the  present  analysis,  incidentally,  is  the  average 
stress-end  shortening  curve  shown  in  Figure  7.  This  solution  corres¬ 
ponds  to  the  deflection  function,  equation  (23),  with  w*  =  0,  A ^  =  0, 
and  a  free  variation  of  the  total  potential  energy  with  respect  to 
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&20>  011(1  n»  Tbe  curve  of  Figure  7  is  not  much  different  from 

that  of  Kempner  in  Figure  4.  The  minimum  equilibrium  postbuckling  stress 

parameter  is  given  by  oR/Et  =  0.195*  This  minimum  value  is  significant 

in  that  it  has  been  obtained  on  the  basis  of  a  buckle  pattern  which  is 

a  first  modification  of  the  classical- theory  chessboard  pattern.  Fran 

,  ,  ft 

the  discussions  contained  in  the  studies  presented  by  von  Karm&n-Tsien, 
Leggett  and  Jones, ^  Michielsen,^  and  Kempner,^  the  implication  is 
always  present  that  the  low  minimum  postbuckling  stresses  obtained  are 
associated  with  the  "diamond"  buckle  pattern  developed  by  von  Barman  and 
Tsien  on  the  basis  of  physical  observation  which  automatically  includes 
the  AQ2  term  in  equation  (23). 

The  difficulty  in  obtaining  an  economically  feasible  solution  beyond 
aR/Et  =  0.334  (on  the  unstable  portion  of  the  curve)  is  attributed  to 
the  fact  that,  for  small  increments  in  the  applied  unit  end  shortening 
(see  Figure  4),  the  buckle  maximum  amplitude,  nondimens ionalized  with 
respect  to  the  wall  thickness  t,  is  growing  very  rapidly  whereas  p. 
and  i]  are  becoming  quite  small  in  magnitude.  The  buckle  pattern  free 
coefficients  4^,  A2Q,  and  AQ2  are  given  in  Figure  8  as  functions 
of  eH/t;  the  free  parameters  ft  and  tj  are  given  in  similar  fashion 
in  Figure  9* 

The  tedious  and,  hence,  cost-consuming  computations  necessary  to  seek  a 
converged  solution  of  the  five  nonlinear  equations  in  A^,  A^,  A q2,  ft, 
and  7]  are  not  justified  because  the  portion  of  the  average  stress -end 
shortening  curve  in  question  represents  application  of  the  von  Karman- 
Dormell  theory  in  a  range  where  it  is  no  longer  valid  for  practical 
considerations;  the  limitation  on  either  the  size  of  the  rotations  or 
the  linear-elastic  constitutive  law  is  being  exceeded, and  the  results 
have  meaning. only  for  shells  of  extremely  high  R/t  ratios.  The 
limitations  of  von  Karman-Donnell  theory  with  respect  to  permissible 
rotations  and  onset  of  inelastic  deformations  for  axially  compressed 
perfect  shells  in  the  R/t  range  of  practical  interest  have  been  pointed 
out  by  Mayers  and  Rehfield."^  In  fact,  in  reference  16  it  is  shown  that 
the  Kempner  solution  of  reference  13  based  on  von  Karman-Donnell  theory 
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is  the  lower  bound  to  a  family  of  R/t -dependent  average  stress-end 
shortening  curves  and  corresponds  to  the  case  of  R/t  -» ■. 

The  present -analysis  solution  shown  in  Figure  4  for  of^Et  >  O.587  does 
not  exceed  the  kinematic  limitations  or.  von  Kano^n-Donnell  theory,  as 
the  maximum  deflections  are  of  the  order  of  the  vail  thickness;  however, 
it  is  of  no  practical  value  because  it  rep:  esents  unstable  equilibrium. 
Above  the  bifurcation  point,  however,  it  is  of  interest  to  note  that  the 
thin-shell  theory  is  attempting  to  produce  thick-shell  behavior 
(axisywmetric  buckling)  because,  as  indicated  in  Figure  9,  T)  -*  0  and 

p  os. 

*  + 

The  von  Karman- Donnell  theory,  Kempner-type,  axially  compressed  perfect 

shell  solution  with  more  terms  Included  in  the  radial  deflection  function 

v  (equation  (23)  with  v*  =  0)  has  been  given  by  Almroth^  (nine 

U  15 

free  parameters)  and  Hoff  et  al.  (15  free  parameters).  Although 
Almroth  succeeded  in  obtaining  a  minimum  postbuckllng  equilibrium  stress 
corresponding  to  oR/Et  =  0-0652,  Hoff  et  al.  succeeded  in  lowering  the 
minimum  postbuckllng  stress  parameter  crR/Et  to  0.0427.  From  the 
overall  results  of  reference  15,  the  authors  are  able  to  conclude  that 
when  the  radial  deflection  is  represented  by  an  infinite  series,  the 
completely  free  parameter  solution  of  the  von  Karman -D onnell  theory  for 
the  axially  compressed  perfect  shell  is  represented  by  v  ->  0,  p  ->  0, 
t)  -» 0,  and  a  minimum  postbuckllng  equilibrium  stress  parameter 
cR/Et  -» 0.  The  buckle  shape  itself  approaches  a  perfect  Yoshimura^ 
developable  surface  -(vanishing  membrane  strain).  Thus,  on  the  basis  of 
the  findings  reported  in  references  15  and  16,  it  took  from  1941  to  1965 
to  establish  that  the  rigorous  stable  solution  of  the  von  Karman-Donnell 
theory  for  an  axially  compressed  perfect  shell  with  (l)  a  finite  number 
of  assumed  buckle -pattern  parameters  absolutely  free  applies  only  to  an 
infinitely  thin  shell  (R/t  -» ® )  and  with  (2)  an  infinite  number  of  assumed 
buckle -pattern  parameters  which  are  absolutely  free  leads  to  the  trivial 
solution  that  the  infinitely  thin  shell  buckles  at  a  vanishingly  small 
stress. 
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Some  final  remarks  are  in  order  relative  to  the  present-analysis  solution 
shown  in  Figure  2  and  that  of  Kempner  shown  in  both  Figures  2  and  4. 
Obviously,  with  decreasing  initial-imperfection  amplitudes,  the  present- 
analysis  curve  of  Figure  2  coalesces  into  the  lower  branch  of  the  Kempner 
curve.  How,  when  the  initial  deviation  pattern  wavelengths  are  larger 
than  the  wavelengths  occurring  in  the  Kempner  solution,  an  increase  of 
the  initial-imperfection  amplitudes  will  not  affect  the  Kempner  solution 
became  the  Kempner  buckle  shape  is  physically  able  to  govern  the  behav¬ 
ior.  On  the  other  hand,  when  the  initial-deviation-pattern  wavelengths 
are  smaller  than  the  wavelengths  occurring  in  the  Kempner  solution,  the 
initial-deviation  pattern  governs  and  grows  in  amplitude;  the  Kempner 
buckle  pattern  is  constrained  from  occurring, which  means  that  an  increase 
in  the  potential  energy  occurs.  Thus,  the  minimum  postbuckling  stress  of 
the  imperfect  shell  is  greater  than  the  minimum  postbuckling  stress 
obtained  from  the  Kempner  solution.  The  implication,  then,  is  that 

stiffening  elements  could  be  placed  on  a  thin  shell  with  spacing/  so 
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arranged  that  the  snap-through  buckle  pattern  cannot  occur.  Thib  behavior 
has  been  illustrated  in  the  laboratory  by  buckling  two  identical  cylin¬ 
ders  which  were  commercially  fabricated  on  an  automated  assembly  line. 

One  cylinder  is  buckled  in  the  unstiffened  configuration  so  as  to  discern 
the  buckle  pattern.  The  remaining  cylinder  is  tested  with  very  light 
stiffeners  located  so  as  to  hinder  the  original  buckle  pattern  from 
forming.  A  significant  increase  in  buckling  load  (almost  50  percent)  has 
resulted  with  a  very  small  weight  penalty.  Photographs  of  the  two 
specimens  tested  to  substantiate  tais  behavior  are. shown  in  Figure  10. 

The  remarkable  decrease  in  buckle  amplitude  of  the  lightly  stiffened 
cylinder  should  be  noted.  The  buckle  pattern  of  the  unstiffened  cylinder 
is  of  interest  because  of  its  similarity  to  the  Yoshimura1^  buckle  pattern. 
The  latter  is  a  polyhedral  surface  of  plane  triangles  developable  from 
the  middle  surface  of  a  circular  cylindrical  shell.  Conventional  cylin¬ 
ders  tested  in  the  laboratory  tend  toward  Yoshimura-type  buckles  only  in 
localized  bands.  The  preformed  plastic  hinges  shown  in  the  test  cylinder 
of  Figure  30  provide  the  mechanism  for  the  development  of  the  Yoshimura 
buckle  pattern  over  almost  the  entire  surface  of  the  unstiffened  specimen. 
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Madsen- Hoff  Solution 


The  behavior  of  postbuckled, axially  compressed,  perfect  circular  cylin¬ 
drical  shells  leading  to  the  trivial  solution  (for  an  infinite  number 
of  terms  in  the  trigonometric  series  representing  the  radial  deflection 

function)  discussed  in  the  preceding  section  has  been  explained  by 
12 

Madsen  and  Hoff.  The  authors  of  reference  12  point  out  that  a  varia¬ 
tion  of  the  total  potential  energy  with  respect  to  all  ouckle -pattern 
parameters  is  unrealistic  in  view  of  the  fact  that  a  cylinder  of  finite 
R/t  ratio  cannot  buckle  into  less  than  2  circumferential  full  waves. 
Thus,  in  reality,  the  trivial  solution  cannot  occur. 

In  the  same  reference,  Madsen  an!  Hoff  then  proceed  to  investigate  the 
behavior  of  practical  shells  in  order  to  estimate  the  magnitude  of  the 
maximum  stress  reached  in  the  presence  of  imperfections.  Although  the 
behavior  of  imperfect  shells  in  axial  compression  had  been  investigated 
by  Donnell  and  Wan,^  Madsen  and  Hoff  discerned  an  error  in  the  Donnell- 
Wan  procedure  and  thus  justified  further  investigation  of  the  problem. 
Fundamentally,  the  criticism  of  the  Donnell -Wan  procedure,  as  noted  by 
Madsen  and  Hoff,  is  that  the  total  potential  energy  should  not  be 
mi  ’imized  vith  respect  to  any  of  the  parameters  describing  the  initial 
deviation  shape.  Some  of  the  results  of  their  analysis  based  on  the 
correct  procedure  are  given  in  Figure  5  for  a  specified  imperfection 
pattern.  The  deflected  shape  is  obtainable  from  equations  (21)  and  (23) 
of  the  present  analysis  when  =  0,  A^q  =  0,  and  AQ2  =  0. 

The  curve  of  Figure  5  for  A^  =  0  is  based  upon  a  fixed  pQ  tiq. 

The  postbuckling  curve  for  A^,  A^,  p,  and  tj  as  free  variables,  as 

discussed  earlier,  has  oeen  obtained  for  the  first  time  in  the  present 

analysis  (see  Figure  7).  The  significance  of  p^  =  1.0  and  t)q  =  0.826 

is  related  to  the  only  reasonable  physical  solution  of  the  linearized 

problem  solved  by  Madsen  and  Hoff.  It  is  of  interest  to  note,  however, 
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that  Loo,  by  making  the  arbitrary  assumption  that  p  =  1,  was  able  to 
calculate  that  7  =  0.826  at  the  bifurcation  point.  Of  further  interest 
is  the  fact  that  the  present  analysis  leads  to  the  values  p  =  1  and 
t]  =  0.826  in  a  free  variation  of  the  nonlinear  formulation  at  the 
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bifurcation  point;  in  addition,  it  can  be  seen  from  Figure  9  that,  on 
the  basis  of  extrapolation,  another  set  of  u  and  values  (approxi¬ 
mately,  p  =  0.5^3  and  J)  =  0.320)  exists  at  the  bifurcation  point . 


The  imperfect  shell  average  stress-end  shortening  obtained  by  Madsen  and 
Hoff  for  =  0.1  is  given  in  Figure  5  as  a  solid  curve.  A  distinct 
maximum  stress  is  evidenced.  This  Maximum,  the  conclusions  of  reference  12 
notwithstanding,  is  not  the  lowest  maximum  stress  to  be  expected.  The 
addition  of  the  term  in  equation  (23)  leads  to  the  dashed  curve  in 

Figure  3  on  the  basis  of  the  present  analysis  when  >.  =  X  .  In 


addition,  for  the  case  when  A£^  =  0.17,  the  inclusion  of  the  A^  term 
removes  the  concept  of  an  initial  maximum  stress.  Thus,  it  is  evident 
that  maximum  stress  calculations  for  reasonably  imperfect  shells  based 
on  the  Jfcdsen-Eoff  procedure  should  include  sufficient  terms  ir.  the 
assumed  deflection  function  to  guarantee  a  converged  solution.  This 
bwumi  that  the  maximum  stress  is  sensitive  not  only  to  the  aiqplitude  of 
the  Initial  deviation  but  also  to  the  waveform  of  the  buckle  pattern. 
Finally,  it  is  important  to  note  that  maximum  stress  calculations  should 
not  overlook  the  possibility  of  yielding  of  the  material.  A  maximum 

strength  analysis  of  axially  compressed  imperfect  shells,  including 

2L 

inelastic  deformation,  has  been  presented  by  Mayers  and  Wesenberg,  '  in 
which  a  modified  version  of  Reissner's  variational  principle  is  utilized. 


CIBCIHAHLY  CURVED  PLATE 

The  analysis  of  the  long  curved  plate  undertaken  in  this  investigation  is 
important  from  the  standpoint  that  conventional  aerospace  vehicle  shell 
structures  are  generally  of  the  semimonocoque  rather  than  pure-monocoque 
type*  That  is,  the  pure  monocoque  thin  shell  discussed  in  the  previous 
section  is  stiffened  by  stringers  and  rings  in  actual  practice  in  order 
to  overcome  the  tendency  of  pure-monocoque  thin  shells  to  buckle  at 
stresses  well  below  the  classical  value  under  direct  compression  and 
bending  loads.  To  date,  despite  extensive  theoretical  and  experimental 
studies,  the  optimum  semimonocoque  shell  for  a  given  loading  (that  is, 
the  sizing  of  stiffeners  and  shell  thickness  and  the  spacing  of 
stiffeners)  cannot  be  designed  on  the  basis  of  theoretical  analysis 
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alone.  It  Is  unusual  that  so  much  attention  and  effort  has  been  devoted 
to  pure-monocoque  shells  vhen  the  fundamental  structural  component 
existing  between  stringers  and  rings  of  semimonocoque  shells  (the  actual 
case)  is  a  curved  plate.  This  is  certainly  true  vhen  discrete  stiffening 
is  utilized  and  conceptually  true  in  the  case  vhen  analysts  consider  that 
sufficiently  closely  spaced  stiffeners  and  rings  may  be  "smeared  out" 
over  the  shell  surface.  However,  even  in  the  latter  case,  information 
is  necessary  to  establish  geometrical  properties  (including  spacing  of 
the  stiffeners),  supposedly  to  preclude  the  occurrence  of  local  buckling 
between  stiffening  elements. 

The  average  stress-end  shortening  curves  appearing  in  Figure  6  for 
Z  *  10  and  20,  respectively,  based  upon  the  radial  deflection  function 
of  equation  (26)  and  the  variational  development  of  Appendix  III  (with 
U  =  1  and  X  =  b),  are  considered  to  be  extremely  significant.  It  is 

y 

to  be  noted  that  vhen  the  plate  is  flat  (R  -*»,  z  ->0),  the  well-known 
behavior  reflects  the  ability  of  the  plate  to  carry  additional  load 
beyond  the  bifurcation  point.  As  shown  in  reference  25,  this  behavior 
is  limited  only  by  the  phenomena  of  eventual  buckle  pattern  change  and 
the  onset  of  inelastic  deformations.  Under  the  assumption  of  elastic 
behavior  in  the  present  analysis,  the  solid-curve  result  in  Figure  6  for 
a  shallow -curved  plate  (z  =10)  suggests  that  catastrophic  snap-through 
buckling  can  be  avoided  and  a  load-carrying  ability  in  the  vicinity  of 
the  classical  value  can  be  maintained.  With  Z  =  20,  the  snap-through  does 
occur,  but  it  cannot  be  considered  catastrophic  from  the  load -carrying- 
ability  standpoint.  Even  when  the  improved  radial  deflection  function 
(equation  27)  is  introduced,  only  about  a  ten  percent  reduction  in  load¬ 
carrying  ability  is  evidenced.  However,  for  Z  =  kO  and  60,  the 
present  analysis  computations  indicate  that  thin-cylinder  behavior 
(catastrophic  snap-through)  governs.  The  lower  bound  to  such  behavior 
for  the  present  analysis  (with  =  0)  is  given  by  the  lowermost  solid 
curve.  This  curve  is  the  same  as  that  appearing  in  Figure  5  for  a 
perfect  cylinder  Wiiich  buckles  and  postbuckles  with  p  =  1  and 
t]  =  0.826. 
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The  primary  question,  of  course,  concerns  the  v eight  penalty  involved  in 
sizing  stiffened  shells  such  that  the  panel  behavior  between  stiffeners 
and  rings  is  given  by  results  similar  to  those  plotted  in  Figure  6. 
Apparently,  in  practice,  cylinders  with  stiffeners  spaced  sufficiently 
close  together  that  the  "smearing  out"  concept  may  be  used  are 
specific Ally  designed  such  that  local  buckling  of  the  sheet  material 

between  stiffeners  Is  prevented.  For  example,  the  recent  preliminary 

2  6  / 

design  procedure  presented  by  Smith  and  Spier  for  stringers  and/or 

ring-stiffened  circular  cylinders  under  axial  compression  commences  with 

the  fundamental  step  that  precludes  buckling  cf  the  curved  sheet  lying 

between  the  stiffening  elements  of  the  overall  cylinder.  As  another 

27 

example,  the  recent  shell  analysis  manual  prepared  by  Baker  et  al. 
recommends  (1)  the  design  formula 
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(28) 


bo  determine  the  buckling  stress  for  the  curved  plate  elements  between 
stringers  and  rings  and  (2)  the  generous  knock-down  values  for  C  extending 
from  40  percent  for  R/t  =  100  to  55  percent  for  R/t  =  500,  the  prac¬ 
tical  R/t-ratio  range  for  the  sheet  material  in  stiffened  cylinders. 

The  coefficient  tj  in  equation  (28)  has  been  introduced  in  reference  27 
as  a  correction  factor  when  inelastic  buckling  occurs. 

To  study  the  implications  of  the  present  results  (see  Figure  6)  as  they 

might  apply  to  achieving  increased  efficiency  of  axially  compressed 
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stiffened  shells,  a  particular  cylinder  tested  by  Card  and  Jones  has 
been  selected.  The  cylinder  geometry  and  elastic  constants  are  given 
in  Figure  11.  It  is  noteworthy  that  even  though  the  effective  R/t 
ratio  of  the  "smeared  out"  cylinder  is  167,  the  "effectively"  thick 
cylinder  buckled  in  a  typical  thin-cylinder  "diamond"  configuration  (see 
Figure  3b  of  reference  28).  This  would  tend  to  indicate  a  significant 
sensitivity  to  initial  deviations  from  a  uniform  circular  shape. 


For  the  cylinder  section  of  Figure  11,  the  Z  of  the  curved  plate 
between  60  stringers  spaced  1.0  inch  apart  is  3«7^*  Now,  when  every  other 
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stringer  of  the  60  stringers  is  removed,  Z  becomes  Ik. 96.  In  an 

attempt  to  provide  for  the  formation  of  nodes  at  each  remaining  stringer, 

29 

the  criterion  given  by  Schildcrout  and  Stein  for  axially  compressed 

curved  plates  vlth  simply  supported  edges  has  been  used  vlth  the  result 

that  the  depth  of  the  stringers  must  be  increased  to  t  =  0.33 6  inch 

s 

The  increased  stringer  material  is  "borrowed"  from  the  30  stringers  that 

have  been  removed.  Next,  the  remaining  material  from  the  removed 

stringers  is  distributed  as  20  rings  of  the  same  geometry  as  the  30 

redesigned  stringers.  With  curved  square  sections,  the  criterion  given 

by  Batdorf  and  Schildcrout^  for  simply  supported  compressed  curved 

rectangular  plates  with  central  chordvise  stiffeners  implies  that  nodes 

should  form  at  each  shell  ring  independent  of  the  stiffness  of  the 

ring  cross  section.  The  ring  spacing  has  been  selected  to  give  u  “  1 

in  every  curved-plate  element.  From  Figure  6,  it  can  be  seen  that  the 

average  stress  remains  in  the  vicinity  of  the  classical  stress  for 

Z  =  I4.96.  Thus,  for  a  6  percent  increase  in  weight  of  the  entire  cylinder, 

the  classical  load  obtained  from  the  buckling  criterion  given  by  Block 

et  al.^1  is  6.55  times  greater  than  the  corresponding  buckling  load  for 

the  original  cylinder  with  60  longitudinal  stringers  and  no  rings.  Even 

though  the  quantitative  results  may  not  be  practically  useful,  the 

qualitative  trend  cannot  be  overlooked.  That  is,  a  judicious  redesign 

with  respect  to  the  rings  should  both  lighten  the  cylinder  and  reduce  the 

buckling  load  to  a  practical  level.  The  attempt  should  be  made  to 

determine,  by  iteration,  the  lightest  cylinder  for  Z  <  20  and  an 

effective  radius -to-thickness  ratio  below  about  250.  For  example,  a 

trade-off  study  could  be  undertaken  with  respect  to  the  effects  of  (l) 

lightening  the  20  rings  and  (2)  doubling  the  spacing  between  the  existing 

rings  and,  if  necessary,  increasing  the  ring  stiffness  through  use  of 

the  material  available  from  the  removed  rings.  In  any  event,  the 

impressive  results  obtained  herein  with  respect  to  the  potential  increase 
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in  efficiency  of  the  cylinder  tested  by  Card  and  Jones  cannot  be  dis¬ 
regarded.  For  eccentrically  stiffened  cylinders,  the  decrease  in 
effective  radius-to-thickness  ratio  minimizes  the  sensitivity  to  initial 
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Imperfections.  The  appearance  of  "diaaood"  buckles  in  the  cylinder  test, 
even  with  the  effective  H/t  ratio  equal  to  167,  suggests  that  ring 
stiffening  can  only  aid  in  alleviating  the  imperfection  sensitivity. 
Finally,  the  sasple  calculations  show  thrt  rings  can  be  added  at  little 
or  no  Increase  in  weight  by  redesigning  the  stringers  vith  respect  to 
section  properties  and  spacing. 


CONCLUDING  REMARKS 


The  extensive  study  presented  has  made  it  possible  to  understand  to  a 
greater  degree  the  long-standing  and  enigmatic  behavior  of  thin  circular 
cylindrical  shells  in  axial  compression.  The  significant  contributions 
of  Donnell, '  von  Karman  and  Tsien,  and  Keener  J  stemming  from  their 
investigations  of  perfect  elastic  shells  have  been  reevaluated  cn  the 
basis  of  the  present  analysis  and  some  new  results  and  interpretations 
given  the  earlier  works.  With  regard  to  the  behavior  of  Imperfect  shells, 
the  analyses  of  Donnell  and  Wan  and  Madsen  and  Hoff12  have  been  studied, 
also,  and  the  results  of  the  latter  authors  have  been  extended. 

On  the  basis  of  the  understanding  of  pure-monocoque  shell  behavior,  an 

investigation  of  the  postbuckling  behavior  of  curved  plates  has  led  to 

significant  results  concerning  the  transition  region  from  flat  plates  to 

highly  curved  plates  and  complete  cylinders.  The  curved  plate  behavior 

is  significant  in  that  the  conventional  thin  shell  utilized  in  aerospace 

structures  design  is  of  the  semimonocoque  (stiffened)  rather  than  the 

pure-mcnocoque  configuration.  Thus,  the  basic  structural  shell  elements 

are  the  curved  plates  existing  between  stringers  and  rings.  In  the 

absence  of  any  known  theoretical  analysis  which  can  predict  the  maximum 

load  in  axial  compression  of  a  given  stiffened  shell,  the  present 

investigation  offers  some  new  concepts  which  appear  to  justify  the 

preliminary  design  of  optimized  conventionally  (stringer-ring)  stiffened 

shells  on  the  basis  of  nonlinear  curved  plate  analysis  in  conjunction 

with  stringer-ring  sizing  and  spacing  based  upon  (l)  the  postbuckling 

results  obtained  herein  for  unstiffened  thin  shells  and  (2)  the  linear- 
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theory  results  of  Schildcrout  and  Stein  and  Batdorf  and  Schildcrout 
for  axially  compressed,  simply  supported,  curved  plates  with  axial  and 
chordwise  central  stiffeners,  respectively.  The  qualitative  results 
obtained  relative  to  achieving  buckling  loads  of  stiffened  cylinders  in 
the  vicinity  of  the  classical  value  cannot  be  overlooked.  Further  re¬ 
search  in  this  area  is  mandatory  because  of  the  desperate  requirement 
for  veig.it  reduction  in  aerospace  structures.  The  latest  design  manuals 
(19 67-68 ;  ’  ’  9  still  contain  cautionary  remarks  to  the  effect,  for 
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example,  that  "A  currently  popular  viewpoint  is  to  consider  classical 
theory  to  be  directly  applicable  to  most  practical  stiffened  shells. 
Nevertheless,  to  account  for  uncertainties  and  to  guard  against  reckless 
extrapolation  into  extreme  parameter  ranges,  it  is  suggested  here  that  a 
knock-down  factor  be  retained  in  the  budding  analysis  of  stiffened 
cylinders.  This  should  result  in  conservative  strength  estimates  which 
can  be  employed  with  confidence  in  the  design  of  actual  hardware." 
Unfortunately,  the  knock-down  factors  suggested  are  very  little  different 
from  those  which  have  been  generously  utilized  in  practice  for  the  past 
25  years  at  the  expense  of  overweight  and  uneconomical  designs.  If 
optimized  design  criteria  are  not  available  for  conventional  shell 
construction,  then  what  hope  is  there  that  such  criteria  will  soon  be 
available  for  shells  in  which  promising  conmcsites  are  to  be  employed 
(for  example,  sandwich,  laminated,  and  fiber-reinforced  materials)? 
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Figure  k. 


Average  Stress-Unit  End  Shortening  Curve  for  Kempner  Radial 
Deflection  Function  (Present  Analysis  Digital  Computer 
Solution  and  Desk  Calculator  Points  Obtained  in  Reference  13). 
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Figure  6.  Average  Stress -Unit  End  Shortening  Curves  for  Flat  and  Curved 
Plates  Covering  the  Range  0  <  Z  <  as. 
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Figure  8.  Variation  of  Radial  Deflection  Pattern  Amplitudes  n , 
ana  A^  With  Unit  End  Shortening  for  Present  Analysis 
Solution  Corresponding  to  Kempner  Deflection  Function. 


Figure  10.  Effect  of  Lightweight  Longitudinal  Stiffeners  in  Preventing 
Catastrophic  Snap-Through  Buckling  Into  a  Yoshimura  Pattern. 
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APPENDIX  I 


DEVELOPMENT  CF  GOVERNING  EQUATIONS  AND  BOUNDARY  CGKSITIOKS 

The  total  potential  energy,  equation  (8),  upon  substitution  of  the 
strains  given  by  equation  (5),  and  subsequent  integration  through  the 
sheU  thickness,  gives  In  expanded  form 


(1-v) 


L/2  2jtR  1  2.1  4  ,1  2  2  1  2  1  4  1  2  2 

J-L/Jo  “5  U'x+Sv'*+5"’*V5  VS 


12  12  12  1 
+  - SK  +7TU,  V,  +U,  V,  V.  +U,  V,  +ttU.  V.  +U,  V,  Wrt  --  U,  V 

2r2  2  'x  'x  'x  'x  0,x  'x  'y  2  'x  *y  'x  >y  0  R  'x 


A 1  3  ±1  212  212  12 

2V,xW0  +2V,y  'x  r^xW,y  2^xW'yV0  -2RW'xW  V,yV'xV0 
x  y  x 


12^  1  1  2 
2  'x  ,y  0,x  'x  'y  0,x  0,y  R  ^x  0,x  2  'y  'y  'y  >y  0,y 

1X3121  12 

-  —  v,  v+~jv,  V-  v-i  v.  w.  v)  -  (l-v)(u.  v,  +iu.  v, 

R  *y  2  ;y  0,  2R  Jy  R  'y  0,  '  'x  'y  2  'x  'y 

y  y 

1  ^12  12^  1 

+  U.  V,  V-  -  =  U,  W+^V.  v,  V.  V+W.  Wn  V,  +=- V,  w~  w . 

*x  *y  0,  R  'x  2  'y  2R  *x  'x  0,  ’y  2  'x  0,  ’y  0, 
y  x  x  y 

1  12121221221  1 

-  —  W,  V-XU.  -TV.  -T-V-.  V,  — p  V.  -7?U,  V,  ~7f  U.  V,  V, 

R  ’x  0,x  4  'y  4  'x  4  0,x  'y  4  'x  0,y  2  'y  'x  2  'y  'x  ’y 

1  1  1  ,  1 

-  7  U|  V.  -oU.  V.  VA  -o  vj  V,  -V-  -o  Vj  WA  V. 

2  *y  0,  *y  2  7y  7x  0j  2  Jx  Jx  7y  2  7x  0,  7y 
x  y  x 

i  nr  ^'/2(i2ItR  2 

-  2  v»x"'x"0,y)1,ixay  +  2  J.L/jo  [(v’=aM''yy) 


-  ^"'^'‘■xx  “'jry-’'’xy))axay 


In  accordance  with  the  potential  energy  principle,  the  first  variation 
of  the  strain  functional,  equation  (29),  with  respect  to  the  displace¬ 
ments  u,  v,  and  v  must  vanish.  Thus, 
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Integration  by  parts  followed  by  appropriate  grouping  of  terms  yields 
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With  the  use  of  the  strain-displacement  relations,  equation  (5),  equation 
(31)  can  be  rewritten  as 
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+  2(l+v)  I  ^xv611  ^  +  - p”  I  (e  +vc V)8v|  dx 

2(1+V}  J-L/2  Xy  0  (l-v2)  J-L/2  y  X  0 


1+9 


( continued ) 


?*R  L/2  L/?  2*R 

+  2Ti~7T  !  7xv6vI  dy  +  D  I  (v,  »-vw,  )&v,  |  dx 

^  -L/2  ‘-L/2  -  **  y0 

r  2xR  L/2  L/2  2nR 

+  D  i  <w'xx+VW'yy)&W'Xl  -  dy  “  D  !  ,  1",  +(2-v)v,  ]&w|  dx 

'0  xx  yy  x  _L/2  --l/2  yyy  ^  0 

P  2nR  .  L/2  ,  L/2 

“  D  I  tw'xxx+^2*v^w'xw^&W  ,  dy  +  - IT  '  t(v.  -4V  )(e°+ve°) 

’  0  ***  ^  -L/2  (1-v2)  °  -L/2  y  °'y  y  x 

PjtR 

+  ^f^(v,+w0  )7°]&h|  dx  +  [  [(v,  +w  )(£x+v€y) 

2  X  0,x  xy  '0  (i-v2 )  *  0  x  °'x  x  y 

/,  \  n  L/2  2itR  L/2 

+  (v  +u  )7  J6w|  dy  f  2D(L-v)v.  5w|  |  =0  (32) 

y  'y  ^  -L/2  xy  0  -L/2 


Nov,  with  the  introduction  of  the  stress -strain  relations  (equations 
(6)),  vritten  for  the  t  idle  surface,  the  variation  of  the  total 
potential  energy  vith  respect  to  the  displacements  u,  v,  and  v  is 
expressed  as 

_  L/2  ,^itR  _  -  n  n 

&  (U)  -  -  f  !  {(r  +N°  )b:  +  (N°  +N°  )8v 

u'v'w  -L/2-0  *x  y'y  y-v  y,x 


y,y  xy,x 


+  N°(v,  +H.  ),  +  2N°  (v,  +wn  ),  +  4  N0  -  DV  v]Bw}dxdy 

y  y  °f'  y  xy  x  °>x  y  R  y 

V 

T  /2 

+  C^x°ySU  +  HySV  +  'y 

2nR  _  2tR  nn  0 

+  2M  lew] |  dx  +  (N  5u  +  N  Bv  +  [N (x,  +w  ) 

xy,x  'o  Jo  y  AX  u,x 

L/2  L/2  2jtR 

+  N  (v,  +w  )  +  M  +  2M  ]&-}  dy  -  J  M  Bv,  |  dx 
xy  y  u,y  x,x  xy,y  _L^2  j-l/2  y  y  0 

2itR  L/2  2nR  L/2 

M Bw,  I  dy  +  2D(l-v)v,  Bw  |  |  =0  (33 

J0  -L/2  y  0  -L/2 
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APPHIBIX  II 


SOLUTIOB  OF  COMPATIBILITY  EQUATIOR 


For  a  given  trigonometric  series  representing  the  initial  and  buckling 
radial  deflect! cm  shapes,  the  procedure  for  integrating  the  compatibility 
equation  (24)  is  as  follows: 

Equations  (21)  and  (23)  can  be  rewritten  as 


F  =  A11  COS  W  C0G  V  +  A20  C0S  2W 


(38) 


and 


^  *  Aqq  +  A^  cos  pax  cos  ay  +  A^q  cos  2pax  +  AQ2  cos  2ay 


+  A^  cos  PgSgX  cos  agy  +  A|Q  cos  2pQa0x 


(39) 


wnere  pQ  =  \y  /\  ,  P  =  %  =  *Ay  /  and  a  *  n/X.y. 

Substitution  o?  equations  (38)  and  (39)  into  equation  (24),  and  perfor¬ 
mance  of  the  indicated  operatic  ;s,  where  v"V*(  )  =  0,  yields 


cos  pax  cos  ay 


F/E  ~  §r  +  (^+1)2  (^2  '  2t  AllA20'2t  A11A02^ 

2  2  \  12  2  2 
"  tAn  cos  2pax  -  Jg  t  p  A^  cos  2ay 


1  /  8tA20  J2. 

+  zrz  v 


32p  '  Ra 

2t  4  AilA20  t  2t  ^  AllA02 

- *  *  1  cos  3pax  cos  ay - * —  »■  ■  ■  cos  pax  cos  3ay 

(9P+1)  (p  +9f 


2  2 
P  t  A„„A, 


'  *20*02  ^0  ftAl 1  2 

5 c°s  2pax  cos  2ay  +  — r - *  — ±=  -2t  (A*  A* 

W  LBaT  U  21 

U5]005  u0a0 


+  aJ1A2Q+A20A 


1  /8tA?0  2  2 

x  cos  a^y  +  - -  I - -  t  A* 

^ v 


(continued) 
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Au)  c“  %a; 


'  2t2 

'x  •  -jr  'Au+2i?iAu>  C0B  2V 


T  2“T7?  ( AiiA2o 'l‘A2cr’lii+AiiA2o  ^  COi:  ^oaox  cos  aoy 

(9ii0+i) 

2  2  2,2 

Moa  aot  C 

rrm  Ao2(A2o+A2o)  cos  2W  cos  2ay 

(Vo4®  > 


t2a2a2(4-n0)2A11(A*L^) 

4[(fia-^0a0)2+(a-a0)^]^ 

t2a2ag(n-40)2A.L1(A*1+A^;i) 

4[(jia+u0a0)2+(a4e0)2]2 


cos  (jia-4_an)x  cos  (a-aA)y 


cos  (ua-Hi0a0)x  cos  (a+aQ)y 


t  a  a0(n+n0)  ^(A^+At^)  s 

— - cos  (na-M^Jx  cos  (a+a^y 

4[(jia-u0a0)  +(a+aQ)  j 
t  a  a  (jx-Hi  )  A  (A*  +A,,; 

— — - -g-.  ■  g  >  -  cos  (na-Hi0a0)x  cos  (a-a  )y 

M(ua+U0a0Ji  +(a-aQ)  j 

_  2  222  .0  > 

2t  8  a0^0All  A20+A20  ,  „  \ 

- s — 5-s -  cos  (toa-2n  an)x  cos  ay 

[(ua-2n0a0)24a2]2  0  0 


2t  a  a0'i0All^A20+A20^ 
[(4a+2nQa0)  +a  ] 


cos  (na+2n  an)x  cos  ay 


0.2  222  / . „  .0  » 

2t  n  a  a0A20(A11+A11) 

— - 3  g3~  =os  (2na-,i0a0)x  cos 

[(2aa-,i0c0)  +aQ] 

2t  ^  a  a0A20(A|1+A11) 

- = — =-* -  cos  (2(ia+u  a n)x  cos  a^r 

[(2^0a0)2+a2]2  ^ 


(continued) 


/ 


— - - T3-  cos  W  cos  (2a-ao,y 


r ,  >2  /r>  N2,2  r0  0  '  0" 

f(^oao}  +(2a"a0)  ] 

2t2u2aca2A  (A*  +A*2  ) 

°  g”  ^  =0.  W  cos  (2a^0)y 
[(40a0)2+(2a+a0)2]2 


Now,  with  the  introduction  of  the  stress  function  F,  equation  (40),  and 
the  initial  and  arbitrary  displacements  vQ  and  w,  equations  (38)  and  (39), 
into  equation  (17),  the  total  potential  energy  is  obtained,  in  terms  of 
the  displacement  amplitude  and  the  buckle  wavelength  parameters,  and  it  is 
expressed  as 


2  U  A 

,crRx  .  ,11 


O  4  2  i 

\2  ^  H  > 


U  =  (— )  +  - ti -  (-— ■  -  nA  A  -nA  A  V  +  —  H  A4 

'Et '  (  2+1)2  k  2  ^V20  ¥11  02j  12b  A11 

„  42. 2  2  4222  ,4222 

1  (  .  _  4  2  .  h  t]  A11Ag0  h  4  AppApg  V  4  A20A02 

2  20  H  H  (9(i2+l)2  (^2+9)2  (^+1)- 


+  12y  (Aii+2AiiAn)  +  ^^2+1^2  ‘■A11'2t10^A11A20+A20A11+A11A20^ 


+  — irr  f8A*  -n  A*  (A*  +2A2  )]2  +  - 2-2 —  (a*  a*  +A*  A0  +A*  A  )2 

12$  1  "20  lo  ir  il  11  ' J  2,,N2  U11  20  20  11  11  02 ' 


,422 

voV 


(w^ 


£  (a  a*  +a  a0  )2  + 
2  n2  v  02  20  02*20 J 


2  2i  \4A2  t  ^.0  \2 
4  40(^^)  All^All+All^ 


(9Pq+1)  ^  xx  xx 

^^>-^)u4(A;i+Aiif 

64[(pv/^-pov/^)2+(/4- 

42420(^0)U4(Aii<)2 


4[(u/4+^0\Au)2+(  /*)+ >Al)C  T  84[(p/^-p Qs/\T+(Jn+  ^4 0) 


'■l\KV  r^gV  XV  1  V  IQ 

)2 


~  ■  v  1  f*Q  v  *o  v  v 

4  2  2  2  ,  A_  O  \2 
^  ^  0A11  A20+A20 
r  /  r"  ^  /”“  \ 2  .  t2 


_  O  O 

v  > 


64[(p/^+b0/T)^+(\/^- /%)  1  [(M-  A-2hn/4n)2+4] 


(continued) 
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4  2  2.£  .  _  n  v2 

40n  Vh  (;':20+a20) 

+  ,,  /—  x2.  i2  + 


4  2  2.2  0  x2 

»  *  T)0A20(A11  All} 


[(^yn+2fi0y^)  +tj]  [(2^7^-^/^)  +Ti0r 


,>  2  2  2  ,  0  v2  4  2  2  2  .  0  x2  4  2  2,2  ,  ...  0  x2 

+  **  71  VWA11+V  +  W  Vo2(A11+A11)  ...  W  Vo2^A11+A1I) 

[(2^i  >/^o)2+ti0]2  [n^c+(2>A+y^')2]‘: 


_2,  2  \2.2 


2/.  4.2 


2/  2 v2.„2 


42.? 


n  (i*  +1)  ^  2l^  ^  ^o+ao2)  W1)  A*n  2Wlo 


48(l-v2) 


3'1-v2) 


48(l-v  ) 


3(l-v2) 


where  tj  =  «2Rt/X,2  =  a2Rt  and  tj-  =  «2Rt/X2  =  a2Rt 

y  u  yo  0 

The  first  term  in  equation  (4l),  that  is,  the  average  stress,  can  be 
expressed  in  terms  of  the  applied  end  shortening  e  and  the  displacement 
parameters  by  utilizing  tne  relation  defining  the  unit  end  shortening 


1  P  L/2 

=  -  i  I  u,  dx 
L  **  -L/2  X 


In  terms  of  the  stress  function  F  and  the  displacements,  equation  (42) 
becomes 

L/2 

■  e  ‘  -  Z  I_L/2(I(F'jy-vI''xX,-(J  '''x«’x',0,>to 

by  virtue  of  equation  (5),  equation  (7)  written  for  the  middle  surface, 
and  equation  (16).  Substitution  of  equations  (38),  (39),  and  (40)  into 
equation  (43)  and  performance  of  the  indicated  operations  yields  (after 
multiplication  of  both  sides  of  equation  (43)  by  R/t) 

•  f  ■  I +  +  4> +  4<tt  +  4 +  +  2aSo4>  w 


The  variation  of  the  total  energy  (4l)  with  respect  to  all  of  the  free 
parameters  yields  a  set  of  nonlinear  algebraic  equations  in  the  variables 
of  the  assumed  displacement,  equation  (39 ),  with  the  prescribed  end 


/ 


shortening  being  the  forcing  function.  The  variation  is  made  indepen¬ 
dently  of  Aqq.  This  parameter  is  evaluated  by  enforcing  the  periodicity 
condition  on  the  displacement  v.  With  the  use  of  equation  ()?).,  equation 
(7)  written  for  the  middle  surface,  and  equation  (16),  the  relation 


51 


APPENDIX  III 


SOLUTION  CF  CURVED-PLATE  PROBLEM 

The  assumed  radial  deflection  shape  in  the  curved  plate  is  assumed  to 
be,  in  general, 

,  rtx  «y  ,  ay  ,  .  .  2t£x  ny 

W  =  tA,,  COS  r—  COS  r4*  +  tAA,  COS  —  tA„-  COS  COS  r4* 

11  A.  A.  01  b  21  a.  A 

x  y  x  y 

However,  in  th*"  present  solution,  the  buckle  aspect  ratio  [i  =  \  /\ 

y  ^ 

is  lssumed  to  be  unity  and  X  =  b/N,  where  N  is  odd  and  L/b  is 


The  corresponding  stress  function  is 

2  2  2  2 

„  „  i  /j.2.2  2Njdc  ,  11  x  2Nny  Z  A21  4N«x 

F  =  E  |-(t  An/32)  cos  -  (-35 - )  cos  cos  — 


2  2  2  2 
t  A21  ..  2N«y  1  A11A21  ...  N*x  t  A11A21  ...  3N*x 

“5 —  cos  b - C - cos  T - 35 - cos  ~b“~ 


t^A  A 

9  .2.  ,  Nrtx  2Nity  11  21  3Nkx  2Nity 

—■  t  A  A  cos  cos  — rf- - -yrri -  COS  — r—  cos  -r-4- 

100  11  d  b  b  0(0  b  b 


N  t  AnA0l  Njtx  (u-i)j 

- 5 - o~o  cos  —  cos  J — r~ 

2[N2+(N-1)2]2  b  b 


N  b  A11A01  Nnx 

- 5 - jr?  cos  ~h~  cos 

2[N2+(N+l)  ] 


Eh  on,  (N+1)*y 

b  cos  b 


2  2 

2N  t  AQ1A21  2N*x  (N”1  ) Jty 

- ■= - 5-0  cos  ~ —  cos  - — r— - 

[4n2+(n-i)2]2  b  b 


2  2  2 

2N  t  A  A  2N?tx  (N+l)ny  .  b  tAll  N*x  Nny 

- = - cos  -r— -  cos  — r-—4-  + - 5-5-  cos  -r—  cos  -rf- 

[4N2+(N+l)2]2  b  b  to2/  b  b 


(continued) 


/ 


bt'xfk  2 

.  *21  2Nnx  Niry  oy  . 

i^7cos~cos^‘%^ 


The  mean  end  shortening  is 


(48) 


1  ?L/2r(ax  °yN  1  2  ,j  a  ,  N^t2  /t2  ,  .2  , 
6  ~  "  L  J.l/2[(T  '  V  E  }  -  2  V'x3dX  =  E  +  -^r~  (AH+4A21} 


(49) 


The  total  potential  energy,  excluding  the  potential  of  the  applied  load 
since  the  end  shortening  is  assumed  constant  for  a  rigid  testing  machine. 


is 


U  =  nVaJ^/64  +  17N  Va^  /128  +  ( 1/64 )( 1/169  +  281/29  JN^^A^  A2 


21' 


11  21 


+  (|f)  +  N^rj2A2  ^  A2.f  02(  1+51 H  )/l6  +  N4n2A^A210?/l6  +  A^/l6 


11  Or" 


IS 


n~oi 


hop  p  p  li  P  P  P  P  P 

K  n  A01A21A  (14-6^)  +  N  +  4A2i/23 


N-i 

+  N^2ai1A01A21^1S/4  “  ^Lw-l)  2  /5*  + 


N-l  N-l 

NnAnA21(-l)  2  ^  +  9Nkr12A^LAoiA2ia5M/8°  -  3Rr)A21A21( -1)^/2^ 


"n"oi"2r  in' 
.2 


-/nA^A^/lw  -  N2T1A21A016/UJt(2N-l)  -  N^A^A^a/h* 


+  N2lA^A0ia/hit(2N+l)  -  8N2T]A01A21A/5n  -  SnSa^/^A^it^N-I) 


01  21 


8n2t!A01A217/5k  +  8N2r]A01A217/5^(2NH) 


N-l 


N-l 

-  2(1+v)[27Nt1A^A21(-1)  2  /800n  -  ^nA^A^ -1)  2  /32n 


(continued) 
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% 


-  ^T^A^aP+e2)/^  -  16N2t)A01A21(72+A2)/2>^. 

-  N2n^1A0132»/l6rt(2N-l)  +  ^nA^A^a2*/^* ( 2N+1 ) 

-  l6N2r)A01A21A2^/25n(2N-l)  + 

2  2 

where  r]  =  jt  Rt/b 

a  =  l/fJJ2 +(N+l)2] 

P  =  i/[n2+(n-i)2] 

/  =  l/tltfP+CN+l)2] 

A  =  1/ (4R2 +(N-l)2] 

♦  -  Id*2  +  Iffl  +  1 

♦  =  4H2  -  4N  +  1 


IN 


■{ 


1  for  N  =  1 


0  for  N  ^  1 


The  boundary  conditions  satisfied  by  the  assumed  deflection  are 


v=0  at  y  =  ±  b/2 
L 

(•  ■L/“ 

cr  dx  =  0 
j-l/2  y 

.L/2 

t  dx  =  0 
-L/2  ^ 
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